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F-Region irregularities; ionospheric; S3-4 Satellite 


The Ionospheric Diagnostic Section of the Upper Air Physics Branch, Sp 
Sciences Division requires understanding similarities and differences be 
tween polar and equatorial spread-F irregularities. This report discuss 
equatorial investigations with major technical, scientific and engineer! 
advances contributing to the current understanding and future plans for 
adaptive communications system design and ionospheric modification and c 
trol. In this report, the S3-4 satellite data analyses is summarized. 
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1.1 


INTRODUCTION 


The Upper Air Physics Branch (Code 4120) at the Naval Research Laboratory 
(NRL) conducts upper air and astrophysical research by means of rockets and 
satellites. This research includes studies of the wavelength distribution of 
solar radiation and the time variations of such emissions as well as the resonant 
scattering of solar radiation by atmospheric species. The Branch also conducts 
a rocket astronomy program to study nighttime ultraviolet and x-ray radiations 
from galactic sources and studies the effects of radiation and particles inci¬ 
dent on the Earth's atmosphere. This data is used to determine the ionic and 
molecular composition, temperature and density of the high atmosphere. The 
Ionospheric Diagnostic Section is involved with measurement of ionospheric 
plasma phenomena; specifically, F-region irregularities at equatorial and high 
latitudes. 

Accumulating information regarding the equatorial F-region ionosphere has 
proven to be important for the communications and ionospheric sciences. This 
interest has brought considerable intensity into equatorial investigations with 
major technical, scientific, and engineering advances contributing to the 
current understanding and future plans for adaptive communications system design 
and ionospheric modification and control. In this report the S3-4 satellite 
data analysis is summarized. 

2.1 EQUATORIAL IONOSPHERIC IRREGULARITIES 

The cause-effect relationship for the equatorial ionospheric irregularities 
was investigated by analyzing the S3-4 and Plumex-1 data. The results of S3-4 
satellite data analyses can be summarized as follows: 

0 Near the nighttime equator there are irregularity structures with 
depletions in plasma density up to 3 orders of magnitude and range in horizontal 
extent from ten's of kms to less than 1 km. These depletions show east-west 
asymmetry with more spectral strength on the western wall of the depletions, 
suggesting more scintillations on the western wall. The mean ion mass, as infer¬ 
red from the ratio square (I /I^^) , is more inside the depletions in most of 
the depletions suggesting more molecular ions inside the depletions. 

° The power spectral analyses of S3-4 data inside the irregular structures 
show that in the intermediate wavelength domain (2km - 25m), the power spectral 
index is 2.4(±0.4). This supports the belief that these irregularities are pro¬ 
duced by collisional Rayleigh-Taylor instability. Simultaneous power spectral 
analyses on electron density and mean-ion-mass have shown a general behavior 
tending to lower power and soften spectra in ion mass fluctuations as compared 
with the fluctuations in total plasma density. 

° The ionospheric irregularities during the equatorial spread-F conditions 
are known to range from ten's (may be hundreds) of kms to 1 meter or even 11 cm. 

To investigate the plasma process(es) responsible for the generation of these 
irregularities which range six orders of magnitude in the wavelength domain, 
the histograms of spectral indices ( obtained by FFT analyses of a number of 
S3-4 orbits) in the different wavelength domains have been drawn. The long 
wavelength domain (50 km >X> 1km) lies near 1.5 while for intermediate 




(lkm >A> 20m) and transitional (100m >A> 10m) wavelengths (A> 10m) two types of 
phenomena have been observed on the power spectra: the peaks at frequencies of 
few hundreds Hz corresponding to wavelengths of few meters and a break in FFT 
near 3 meter wavelength. These are explained on the basis of excitation of 
drift wave and lower-hybrid instability respectively. 

The total spectra ranging from ten's of kms to few meters is explained by 
the concept of cascading of irregularities or heirachy of instabilities. The 
equatorial spread-F starts at wavelengths range of hundreds of kms to tens of km 
by either TID s or gravity waves. The longer wavelengths generate the condition 
which are ripe for the onset of Rayleigh-Taylor (R-T) instability. The gradient 
scale sizes for R-T instability is near 20 km. Once R-T instability has taken 
over, it produces the gradient scale sizes of few hundred meters which start 
gradient driven drift instability in the transitional wavelength domain. This 
phenomena of cascading of irregularities to shorter and shorter sizes produces 
conditions for excitation of drift waves, parametric instabilities and lower- 
hybrid instabilities. 

The most significant results of the investigations have been the two- 
dimensional picture of depletions and the understanding of different plasma 
processes in the different wavelength domains of ionospheric equatorial irregu¬ 
larities. 

3.1 LONGITUDINAL VARIATIONS IN EQUATORIAL IRREGULARITIES 

One of the advantages of STP/S3-4 was that it was sun-synchronous and 
polar orbit satellite which made it possible to study the longitudinal variation 
in equatorial irregularities, as well as global coverage. During the six month 
(March-September 1978) life of the satellite, about 24 percent orbits of the 
total number of orbits which covered the nighttime equator observed the ionos¬ 
pheric irregularities. Although the irregularities can be present at almost 
any longitude, 50 percent of the observed irregularities lie in the narrow 
longitude zones of 0-15° and 165-195°E. 

4.1 CALCULATION DEVELOPMENT 

To perform the above tasks, a large number of routines were developed for 
both PDP 11/10 and PDP 11/34 minicomputers. A partial list of these routines is 
given below: 

° To read the data from tapes provided by Aerojet for S3-4 data. 

° Calculation and the plot routines for the saturation currents 
(I (E), X_(I)) and the ratio square [(I R (E)/I (I)) 2 ] for S3-4, 
Plumex-1, and DNA/PLACES data in various formats. 

0 Plot routines for I-V characteristics for S3-4 data and the temperature 
calculation routine. 

* Calculation of gradient scale length both for S3-4 and Plumex-1 
data and subsequent plot routines. 

° FFT analyses routine along with the addition of window and filter. 

8 Calculation and plot routine for a chapman layer used in BIME 
simulations and SAR applications. 
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Satellite and rocket 

observations of equatorial spread-F irregularities: 
a two-dimensional model 


E. P. Szuszczewicz. M. Singh* and J. C. Holmes 
E. O Hulburt Center for Space Research. Naval Research Laboratory, Washington, DC 20375, U.S.A. 

Abstract—Recent rocket and satellite measurements of equatorial F-region irregularities have been 
able to resolve wavelengths comparable to the meter-size sensitivities of the Jicamarca and Allair 
radar backscatter techniques. In a July 1979 rocket campaign at the Kwajalein Atoll, vertical profile 
measurements by ‘in situ' plasma probes showed the F-region marked by a number of large scale 
plasma depletions, each having its own distribution of smaller scale irregularities and a trend toward a 
co-location of the more intense irregularities with positive gradients of larger scale features. Similar 
measurements on the S3-4 Ionospheric Irregularities Satellite have shown large scale depletions (1-3 
orders of magnitude) with east-west asymmetries that point toward the western wall as the sight for 
the more intense plasma density fluctuations. The combined rocket and satellite measurements provide 
a two-dimensional model of macroscopic F-region depletions with small structures tending to develop 
more readily on the top and western boundaries. The model and associated power spectral analyses is 
in concert with a developing catalog of radar observations and the predictions of numerical simulations 
which employ the Rayleigh-Taylor instability as the primary mechanism for the generation of 
intermediate wavelength irregularities. 


I. INTRODUCTION 

Recently, there have been considerable advances in 
the understanding of equatorial spread-F as a re¬ 
sult of improved ground based radar observations 
and coordinated in situ measurements. The long 
standing Jicamarca radar results [e.g. Woodman 
and LaHoz (1976) and associated references] have 
been studied in greater detail by the expanded 
capabilities of the Altair radar facility in the Kwa¬ 
jalein Atoll (Tsunoda and Towle, 1980). These 
more recent results have investigated plume de¬ 
velopment in an east-west cross-section and pro¬ 
vided some tentative identification of plume ele¬ 
ments with local F-region plasma depletions 
(Tsunoda, 1980a, b). 

Early efforts to examine the exact relationship 
between radar plumes and ionospheric depletions 
by performing simultaneous in situ and ground- 
based radar observations (Kelley el al., 1976: 
Morse et al., 1977) were limited to conditions of 
bottomside spread-F and required extrapolations in 
space and lime to establish correlations. These 
frustrations, as well as two unsuccessful attempts at 
Kwajalein in 1977 and 1978, were finally relieved 
by the successful DNA/PLUMEX campaign in 
1979 (Szuszczewicz et al., 1980a) which defini¬ 
tively established the co-location of intense radar 
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returns with the upper region of a topside F-layer 
depletion. The PLUMEX results also provided the 
first vertical profile of large-scale plasma depletions 
with superimposed distributions of smaller scale 
irregularities. 

In an attempt to establish a 2-dimensional east- 
west in situ profile that complements the develop¬ 
ing catalog of Altair observations, we present in 
this paper a composite of rocket and satellite obser¬ 
vations, from the PLUMEX campaign on the one 
hand and the S3-4 satellite on the other. These 
results are summarized in subsequent sections and 
synthesized in a way that suggests a two- 
dimensional model of macroscopic F-region deple¬ 
tions with smaller scale structures tending to de¬ 
velop more readily on the top and western bound¬ 
aries. 

Z. ROCKET-BORNE MEASUREMENTS 

The coordinated measurements of equatorial 
spread-F conducted during July 1979 at the Kwa¬ 
jalein Atoll involved the launch of two in¬ 
strumented rocket payloads designed to probe the 
detailed in situ structure of the turbulent ionos¬ 
pheric plasma. The first launch operation 
(PLUMEX I , 17 July 1979; 0031:3025 LT) was 
conducted during the late phase in the development 
and decay of spread-F. The associated results and 
discussions of ground-based and rocket-borne diag¬ 
nostics have been presented by Szuszczewicz et 
al. (1980a, b). 
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1purposes of comparison with subs. .juv m ptc 
°f S3-4 satellite data we show lieu- in 
i.,,. 1 the uplejt measurements of rclativi ds rsttj u» 
presented by correlated ion- (/„') and electron 
saturation /.'(•V ) currents collected bv a pair of 
0IV board pulsed-plasma-probes (Hoi mi s and 
S/rszcztwicz, 1975: Szuszczewicz and 
Hoi Mts, 19801 in the PLUMEX 1 operation. The 
ordinals has a linear scale for time-after-launch 
ssitli altitude superimposed at 50 s increments. (Be¬ 
cause ion and electron saturation currents have 
sunificanily different sensitivities to velocity, sheath 
and magnetic field effects, variations in /„' and 
/ r '(V*) not mutually corroborated were attributed 
to the various aspect sensitivities and excluded 
from Fig. 1. This approach facilitated analysts, re¬ 
duced computer time, and established credibility in 
ihc interpretation of the curves as relative electron 
density profiles.) 

The results in Fig. 1 show that a number of 
major depletions (AN./N,°*;0.9) were distributed 
throughout the F-region. Each of the large scale 
depletions (identified alphabetically) has its own 
distribution of irregularities, illustrated in Fig. 2 by 
the expanded view of regions C, D, H and /. It is 
clear that C is not a single narrow biteout but a 
collection of rather large irregular structures ex¬ 
tending over a total altitude domain of about 
12 km. (Vehicle velocity in region C was 2.4 km/s.) 
To develop a quantitative view of irregularity fluc¬ 
tuations observed in the F-region, contiguous 
linear detrends over 4-s intervals were executed 
throughout the entire upleg trajectory. The varia¬ 
tions about those linear detrends were then plotted 
in Fig. 1 as ‘Irregularity Intensity,’ with a maximum 
relative scale of ±4. A fluctuation as great as ±4 
approximately represents a ±80% fluctuation about 
the linear detrend. (Correlation of these results 
with macroscale gradients and Altair backscatter 
contours are discussed by Szuszczewicz et ai, 
1980a). The results show that the most intense 
irregularities occurred on the bottomside gradient 
(region C) with corresponding measurements at all 
other altitudes at a much lower level. We note that 
the fluctuations in the largest depletion (region 
H-l) are smaller than those at C. Furthermore, the 
fluctuations at C, D, E and I are more intense than 
at adjacent locations; and C,D,E and I are co¬ 
located with positive density gradients of large scale 
features. 

3. S3-4 SATELLITE RESULTS 
While rocket-borne instrumentation can provide 
vertical profiles of irregularities, a circular or near¬ 
circular orbiting satellite with high resolution in- 


ittiuiiiv.nl.tlniM is required tv assemble irregularity 
plane. The S3-4 satellite tunics just such an experi- 

The S3-4 ionospheric irregularities experiment 
employs a pair of pulsed-plasma-probes (Holmes 
and Szuszczewicz. 1975: Szuszczewicz and 
Holmes, 1980) on a polar (96.4° inclination), sun- 
synchronous (2230 h LT, equatorial crossing), F- 
region orbiting satellite. The experiment provides 
direct measurements of the ionospheric state 
((V„ Tits condition of irregularity (6N,), and 
associated electron density fluctuation power 
spectra [P(fc)] with 5-20 m resolution. 

Figure 3 presents a sample of relative electron 
density data collected during a nighttime equatorial 
crossing with the satellite orbiting at an altitude of 
240 km with a velocity of 7.53 km/s. The velocity 
component pcrpendici''ar to the geomagnetic field 
is 2.0 km/s. 50 s of data are displayed covering a 
total horizontal extent equal to 366 km, moving 
from east-to-west with increasing time. 

While the data represents a transit from the 
southeast to the northwest, our discussions and 
associated interpretations will assume that all varia¬ 
tions along the magnetic field are negligible. We 
therefore interpret our results in the framework of 
a two-dimensional system, perpendicular to the 
magnetic field. This assumption is similarly em¬ 
ployed in theoretical analyses of equatorial spread- 
F (e.g. Scannapieco and Ossakow 1976) and 
supported by 6300 A airglow intensity measure¬ 
ments of F-region depletions which show signific¬ 
antly greater variation in east-west extent than 
north-south (Weber et ai, 1978). Our discussions 
will also assume temporal variations to be nonexis¬ 
tent during the satellite traversal across any given 
depletion (topically 10 s or less in Fig. 3). We 
consider this a reasonable position since saturation 
times for unstable Rayleigh-Taylor modes are of 
the order of 4000 s (Keskinen et ai, 1980a) and 
F-layer vertical drift velocities are generally in the 
10-20 m/s range (e.g. Szuszczewicz et ai, 1980). 
In neither case should our interpretations be in¬ 
fluenced by a 10 s time lapse as the satellite passes 
from the eastern to the western boundary of a 
depletion. 

The data reveal four large scale depletions rang¬ 
ing up to 98 km wide and a factor of 300 in depth, 
with smaller scale structure visible down to a frac¬ 
tion of a kilometer. The upper portion of Fig. 3 
presents ‘irregularity intensity' as measured by rela¬ 
tive fluctuations about continuous 0.27 s (2.1km) 
linear detrends. 

Further discussion is facilitated by identifying 










Fig. 3. The electron density variations as measured by the electron probe aboard S3-4 on rev. No. 

2123. On the top is the irregularity intensity (5//I) determined over contiguous 2.1km intervals 
throughout the depletions. 

certain features in Fig. 3. First, there are clearly further explored through power spectral density 
defined regions of undisturbed background ionos- analyses. We present in Fig. 4 just such results for 
phere, marked alphabetically A through D, the each of the boundaries in Fig. 3 (IE and 1W refer 
smoothness of the relative density and the corres- respectively to the eastern and western boundary of 
ponding 0% fluctuations attest to their undisturbed depletion number 1). Power spectral analyses are 
nature. The eastern boundary is defined as the presented across the boundaries of each of the four 
region of density gradient moving westward from depletions with spectral indices (»]■ in a spectral fit 
an undisturbed domain to the least lower bound of to P = P u /“") ranging from 1.9 to 2.3. (These indi 
relative plasma density in the depletion. Everything ces are in agreement with the work of Kfskinkn et 
to the west of that minimum is defined as the al. (1980) where numerical simulations predict n 
western boundary. Admittedly these definitions 2.0 to 2.S for horizontal irregularity structures per- 
leave no room for a region that might be called the pendicular to B. This result runs parallel to an 
‘depletion center,’ but depletions 3 and 4 suggest earlier comparison (Szt- s/c zt wicz and Hoi.mi s. 
that there is no easily defined ‘depletion center.' 1980; Kesmnen el al. L980b| which showed that 
This is also borne out by other sets of P J data. the power spectral index across region C in 
Focussing on depletions 3 and 4, we see that the PLUMEX l also agreed with the predictions of 
irregularity intensities are 2 to 3 times larger on the Keskinen el al. (1980a) in the vertical plane 1 
western boundary than on its eastern counterpart. More important to the discussion of east-west 
This same relationship is true in depletions 1 and 2, asymmetries is the ratio of the spectral strengths 
but only after a qualification that suggests that 1 (P 0 ) across the two boundaries. By defining Pi as 

and 2 are halves of a larger depletion bounded by the P 0 (West)/P 0 (East), we find the ratio Pi to ex- 
A and B. This is supported in part by the nonexis- tend from 1.4 to 11.3, i.e. the irregularity spectral 
fence of a quiescient ionosphere between the two. strength on the western wall is 1.4 to 11.3 times 
When viewed from this perspective the western more intense than on the eastern counterpart, 
boundary is approximately twice as intense in ir- These results are in concert with the radar obser- 
regularitv intensity as the eastern boundary vations of Tsunoda (1979) which show that the 

Signatures in irregularity strengths and relation- bottomside backscatter strength is often asymmet- 
sliips to plasma instability mechanisms c m be ric in the east-west plane, with the western side of 









Fig. 4. Power spectral analyses across the boundaries of the four depletions (1-4) in Fig. 3. The 
number in the upper portion of each | tel is the spectral index N. in a P = " fit to the data. The 

value P, shown in each right hand pan orresponds to the ratio of spectral strengths P, ■* P 0 (western 
boundary) * P 0 (eastern boundary). 



the plumes being the more intense. The combined 
observations support a model of E-W asymmetry 
which allows for a neutral-wind-driven instability- 
crowth-rate enhancement on the western side of a 
rising bottomside F-region depletion (Tsunoda, 
1979). 

4. CONCLUSIONS 

The combination of recent rocket and satellite 
data allows for the development of a two- 
dimensional empirical model of equatorial deple¬ 
tions and associated irregularities. Elements in the 
data have led to the following conclusions: 


(1) Equatorial depletions are macrostructures 
within which much smaller scale irregularities are 
imbedded; 

(2) These smaller scale irregularities (in the ver¬ 
tical plane) tend to derive their energies from posi¬ 
tive density gradients on the topside of local deple¬ 
tions; 

(3) Vertical and horizontal power spectral anal¬ 
yses display indices from 1.8 to 2.5, a result which 
is consistent with Rayleigh-Taylor instability gener¬ 
ation of intermediate wavelength irregularities dur¬ 
ing the occurrence of equatorial spread-F (Keski- 
men er al. 1980a); 
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(4) When viewed horizontally, bottomside deple¬ 
tions have an east-west asymmetry with the more 
intense fluctuations and spectral strengths being 
observed on the western boundary. This result is 
consistent with the radar measurements of 
Tsunoda (1979) and scintillation observations of 
Livingston et at. (1980): 


(5) The combination of all result:, leads to a 
two-di. nsional model of macroscopic F-reeion 
depletions with smaller scale structures ties eloping 
on the top and western boundaries. This model and 
associated spectral indices fully supports the com¬ 
putational work of Keskinen et at. (1980a). 


REFERENCES 


Holmes J. C. and Szuszczewicz E. P. 

Kelley M. C., Haerenoal G., Rattler H.. 
Valenzuela A., Balsley B. B.. 

Carter D. A., Ecklund W. L.. 

Carson C. W„ Hausler B. 
and Torbert R. 

Keskinen M. J., Ossakow S. L. and 
Chaturveoi P. K. 

Keskinen M. J., Szuszczewicz E. P., 

Ossakow S. L. and Holmes J. C. 

Morse F. A., Edgar B. C., Koons H. C., 

Rice C. J., Heiktula W. J„ 

Hoffman J. H.. Tinsley B. A., 

WtNNINGHAM 1. D.. CHRISTENSEN A. B., 
Woodman R. F., Ponalaza J. and 
Teizeira N. R. 

Scannapieco A. J. and Ossakow S. L. 
Szuszczewicz E. P„ Tsunoda R. T., Narcisi R. 
and Holmes J. C. 

Weber E. J.. Buchau J.. Eathf.r R. H. 


1975 Rev scienl. Instrum. 46, 592. 

1976 Geophys. Res. Lett. 3, 448. 


1980a J. geophys. Res. 85, 1775. 
1980b J. geophys. Res. Submitted. 
1977 J. geophys. Res. 82, 578. 


1976 Geophys. Res. Lett. 3, 451. 
1980 Geophys. Res Lett. 7, 537. 

1978 J. geophys. Res. 83, 712. 


Woodman R. E. and La Hoz C. 


1976 J. geophys. Res. 81, 5447. 


Reference is also made to the following unpublished material: 


Livingston R. C. Weber E. and Bauchau T 

Szuszczewicz E. P. and Holmes J. C. 
Tsunooa R. T. 


1980 Paper presented at International Symposium on 
Equatorial Aeronomy, Aguidilla, P.R. 

1980 NRL Memorandum Report 4289. 

1979 Topical Report 4, SRI. International. 


i 





.. 


; j 

































associated with a density gradient and/or current-driven plasma instability. We have 
made preliminary analyses of the density fluctuations and FFT power spectra for 
evidence of characteristic scale sizes and power law dependence. Scale sizes from 
10-300 km are clearly evident in the irregularities; power law fits to the spectra have 
spectral indices in the range -1.5 to -1.7. For the diffuse aurora our results suggest 
support for a recent numerical study of the nonlinear evolution of the current 
convective instability. 
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THE STP/S3-4 SATELLITE EXPERIMENT: 

HIGH LATITUDE LARGE SCALE DENSITY IRREGULARITIES 


INTRODUCTION 

The polar regions of the earth's ionosphere are known 
to exhibit a wide range of plasma kinetic interactions that 
are very dynamic and complex. The plasma kinetic processes 
have interaction scale sizes that range from the order of 
auroral longitudinal dimensions down to particle gyroradii. 
Although the phenomenology and morphology of auroral processes 
are well known, the physical mechanisms are not well under¬ 
stood. Plasma density irregularities in the auroral F- 
region have been detected through ground-based scintillation 
observations (Aarons, 1973; Fremouw et al., 1977; Buchau et 
al., 1978; Rino et al., 1978), satellite measurements (McClure 
and Hanson, 1973; Dyson et al., 1974; Sagalyn et al., 1974; 
Phelps and Sagalyn, 1976; Clark and Raitt, 1976; Weber and 
Buchau, 1981), and rocket observations (Kelley et al., 

1980). Several physical mechanisms have been suggested to 
explain the high latitude irregularities. A recent review 
(Fejer and Kelley, 1980) lists three general sources of the 
irregularities: 1) particle production (especially by low 
energy precipitating electrons), 2) electrostatic turbulence, 
and 3) plasma instabilities. Except in extremely simple 
situations, it is likely that all three mechanisms will 
contribute to the spectrum of irregularities. 

In this report we will examine two cases of high latitude 
measurements of electron density irregularities obtained in 
the nightside auroral F-region with the pulsed plasma probe 
experiment on the STP/S3-4 satellite. These two cases 
suggest the importance of plasma instabilities associated 
with particle precipitation, i.e, field-aligned currents, 
and illustrate the variety of density irregularities en¬ 
countered at high latitudes. The first observation we 
discuss is that associated with a quiescent aurora, in which 
the structure of density irregularities is relatively well 
defined on several scale lengths, from 10 km to about 300 
km. These scale lengths are significant for some recent 
theoretical studies on the applicability of the current 
convective instability in diffuse auroras. The second 
observation we discuss was obtained during a more dynamic 
state of the aurora and serves to illustrate some of the 
complexity associated with possible gradient- and/or current- 
driven plasma instabilities in discrete auroras. 

MEASUREMENTS AND ANALYSIS 
Plasma Density Measurements 

The polar-orbiting STP/S3-4 satellite is instrumented 
with a pair of pulsed-plasma-probes that allow high resolution 
measurements of electron density N g , temperature T g , density 
fluctuation power spectra P N (k), plasma potential 6 and 
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variations in the mean ion mass 6M./M. [Szuszczewicz et al., 
1981]. We briefly review the experimental operation only 
for the particular data mode of interest here. In the low- 
data-rate mode, the double probe system sampled electron- 
and ion- saturation currents simultaneously at 100 samples 
s 1 . Under conditions of fixed-payload potential and stable 
velocity-field and magnetic-field aspect (see e.g., Szuszczewicz 
and Takacs, 1979), the density N and density fluctuation 
AN are directly proportional to e the probe saturation current 
I and its fluctuations AI . The STP/S3-4 satellite was in 
a e sun synchronous orbit with an inclination of 96.5° 
and crossed the equatorial plane at local times of 2230 and 
1030. The altitude ranged from 160 to 260 km. In the 
discussion below of irregularity scale lengths, we will use 
an average spacecraft velocity of 8 km s 1 to convert time 
variations to spatial variations, using AL=vAt, where AL is 
the spatial length, v is the spacecraft velocity, and At is 
the sampled time interval. 

Quiescent Aurora 

About five minutes of saturation currents, I and I., 
obtained on orbit 390 are shown in Figure 1. These measurements 
were obtained in the night side auroral F-region over the 
south pole at an altitude of 260 km with the spacecraft 
moving from the south polar cap region toward the nightside 
equator (left to right in Figure 1). The probe currents 
clearly show that irregularities were detected on both 
electron and ion probes. The detailed tracking of one probe 
by the other indicates that the probe system was working 
properly and that the irregularities were due primarily to 
variations in the ambient electron and ion densities and not 
the result of secondary effects such as changing spacecraft 
potential, aspect sensitivity, etc. In this preliminary 
report, only the electron density irregularities will be 
discussed, and we simply note the correlative behavior of 
the ion density irregularities. 

The large depletion in electron density (almost an 
order of magnitude) centered at 1840 UT corresponds to a 
diffuse aurora. The geomagnetic Kp index for the period of 
this measurement was 1+, indicating relatively quiescent 
auroral conditions. The ion density shows a similar depletion. 
Both the electron and ion depletions extend over a time 
interval of about 25 s, corresponding to a distance of about 
200 km along the approximately south-north direction of 
spacecraft motion. This aurora was also detected by Huffman 
et al. (1980) with the vacuum ultraviolet (VUV) experiment 
on board the same spacecraft. A detailed comparison of the 
density measurements in Figure 1 with the VUV measurements 















(Fig. 15 in Huffman et al.) indicates that a good correlation 
exists. A few discrete arcs detected by the VUV measurement 
in the time interval 1839-1840 UT appear to be correlated < 

with some of the density structure. Equatorward (to the 
right) of the large depletion, the electron density shows a 
relative enhancement with a scale length of about 300 km. 

A more detailed plot of the electron density depletion 
of Figure 1 and its contiguous domains is shown in Figure 2. 

(The expanded interval is indicated by the arrows above the 
electron current trace in Figure 1.) In this plot every 
data point sampled at 100 s * is shown and the density scale 
on the left of the figure is approximate. Within the depletion 
at least three well-defined density variations with nearly 
periodic behavior are indicated. Assuming the auroral structure 
to be relatively stable, the width of the three density 
structures along the spacecraft trajectory is about 20 km 
and the large scale depletion has a width of about 200 km. 

On either side of the large scale depletion the density 
irregularities generally appear to have high and low frequency 
components with the low frequency components having scale 
widths comparable to the widths of the density structures 
within the depletion, i.e., about 20 km. 

The Fast Fourier Transform (FFT) power spectrum of the 
density irregularities in each of three intervals are shown 
in the insets to Figure 2. The data time interval over 
which the FFT was calculated is indicated below each spectrum. 

The Nyquist frequency for each spectrum is 6.25 Hz, corresponding 
to a minimum scale size of about 0.6 km. Each spectrum is 
normalized to the maximum power calculated from the FFT. 

The straight line drawn through each spectrum is the least- 
squares power law P=P f n , where P/P is the normalized 
power, f is the frequency, and n is°the spectral index. 

The spectra are calculated without smoothing, except for 
the removal of the lowest frequency (DC) component by applying 
a linear detrend. 

The first power spectrum (on the left) suggests a k- 
dependence of k 1 * ** 8 , where k = 2ir/X is the wave number and 
X is the wavelength of the spatial irregularity. The value 
of the index is consistent with previous reports (Dyson et 
al., 1974; Sagalyn et al., 1974; Phelps and Sagalyn, 1976; 

Kelley et al., 1980) on high latitude irregularities where 
scale sizes from less than 1 km up to 200 km have been 
considered. In the second power spectrum, taken within the 
large scale depletion, the spectral index implies a k 1,58 
Variation in spatial irregularity. The low frequency 
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Fig. 2 - An expanded time-scale plot of electron density irregularities in 
the diffuse aurora of Figure 1. The power spectra indicate the frequency 
dependence of the least-squares power law. Frequency f is directly propor¬ 
tional to wave-number k, assuming the irregularities are spatial structures. 









spectral range 0.1-0.3 Hz corresponds to the well-defined 
structures of about 20 km half-width in the density-vs- 
time plot. The spectral power in the 20-km structures 
inside the large density depletion is about three times 
greater than the power in similar scale size structures 
outside the depletion, a fact not obvious in the power 
spectra because of the normalization. The third power 
spectrum shows a power law dependence with a spectral index 
of -1.69. 

An examination of the density-vs-time plot of Figure 2 
shows that approximately 20-km scale size structures are 
clearly evident both inside and outside the large depletion. 

A possible source of these 20-km density irregularities is 
suggested by recent theoretical studies on the current 
convective instability in diffuse auroras (Ossakow and 
Chaturvedi, 1979; Keskinen et al., 1980; Huba and Ossakow, 
1980). The current convective instability is proposed as a 
mechanism for generating large scale (<\/ 10-50 km) density 
irregularities in the diffuse auroral F-region. The current 
convective instability may thus provide an explanation for 
the sheet-like density irregularities that have been detected 
by Fremouw et al. (1977) and Rino et al. (1978) in their 
analysis of enhanced scintillations at auroral latitudes. In 
a numerical analysis, Keskinen et al., (1980) considered the 
nonlinear evolution of the current convective instability 
with an initial density gradient scale size of 50 km. They 
obtained a spectrum of irregularities with a spectral index 
of -1.7 in the north-south direction. This scale size and 
spectral index are similar to the scale lengths and to the 
spectral indices shown in Figure 2. Thus we suggest that 
the observations provide experimental support for the current 
convective instability as a possible mechanism for the 
density irregularities associated with the diffuse aurora. 

Dynamic Aurora 

The plasma probe measurements of Figure 3 illustrate 
the case of a dynamic aurora observed on orbit 244 in the 
southern night-side auroral F-region. The K index for this 
period was 3+, indicating that moderately active geomagnetic 
perturbations were occurring. As in Figure 1, the electron 
and ion saturation currents were sampled in the low- 
data-rate mode (100 s~ l ). The two currents track each other 
very well so that we are confident of proper measurements 
directly proportional to the ambient plasma density. The 
outstanding features of the density measurements in Figure 3 
are the sharp transition in density occurring at about 
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1915:15 UT and the large peak-to-peak variations between 
1916 and 1917 UT. The sharp transition in density is more 
apparent in the electron current than in the ion current 
because of noise interference occurring in the ion probe 
before about 1915:10 UT. A detailed comparison of the two 
currents shows that the transition itself is unaffected by 
the noise and the two currents track each other almost 
exactly within the transition. 

A detailed plot of the electron current in the time 
interval indicated by the arrows in Figure 3 is shown in 
Figure 4 with an approximate density scale to the left. The 
sharp transition begins at about 1915:17 UT and consists of 
large irregular oscillations over the transition width of 
about 1.5 s. The density jump associated with the transition 
is about a factor of 4. Immediately behind the transition 
there are at least four distinct oscillations with an average 
period of about 2 s. These oscillations seem to decay 
exponentially, as indicated by the dashed line (on a log 
scale, an exponential curve will appear as a straight line). 

The VUV observations on the same spacecraft indicate that 
auroral emissions were detected at the same time the density 
transition occurred (R.E. Huffman, private communication) 
suggesting that the density jump is associated with a discrete 
auroral arc. The observation of a sharp transition suggests 
that a moving boundary has been detected and we conclude 
that the density jump is associated with the surge of an 
auroral arc. The surge direction cannot actually be determined 
from our observations; generally the movement is westward 
and poleward (Akasofu, 1968). The sharp transition in 
density implies that gradient-driven plasma instabilities 
are likely to occur; with our interpretation of a surging 
auroral arc, the decaying oscillations behind the density 
jump would represent spatial oscillations left behind by the 
movement of the source of instability, the surge front. The 
electron density increase by a factor of 4 across the surge 
front probably results from the enhanced ionization due to 
precipitating particles. Similar ionization enhancements of 
the F-layer associated with discrete auroral arcs have been 
observed with ionospheric sounders (Pike et al., 1977). The 
surge front very likely has field-aligned currents associated 
with it (Iijima and Potemra, 1978), so that current-driven 
plasma instabilities are also likely to be present and may 
contribute to the density oscillations. It is possible that 
the large and irregular variations seen at later times 
(after 1915:55 UT) also may have had their origin in the 
surge of the auroral arc. If the irregular density structures 
behind the surge front are quasi-stationary, then a wide 
range of scale lengths is evident, from about 4 km to 70 km 
and longer. In contrast, the region ahead of the surge 
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Fig- 4 - An expanded time-scale plot of the electron density for the 
aurora of Figure 3. The sharp density jump is identified as the surge 
of a discrete auroral arc. An exponential decay for the well-defined 
oscillations behind the surge front is suggested by the dashed line. 














front has a distinctly different appearance, with only low 
amplitude, high frequency variations. This region ahead of 
the surge front may be a relatively quiescent ionosphere 
about to be perturbed by the surging auroral arc. The 
different character of the irregularities on either side of 
the density jump that we have identified as a surging 
auroral arc is additional evidence that the density jump is 
a dynamic transition between two states of the F-region 
ionosphere. 

The power spectra for three time intervals are also 
shown in Figure 4. Just ahead of the density jump the 
spectrum of the irregularities shows that the slope of the 
least-squares power law is relatively flat, suggesting that 
thermal fluctuations dominate the spectrum and the ionosphere 
is in a quiescent state. The spectrum immediately behind 
the density jump shows a power law with steeper slope, 
indicating the presence of enhanced power at low frequencies. 
The well-defined decaying oscillations behind the density 
jump appear in the spectrum as the component centered on a 
frequency of about 0.3 Hz. Under the assumption that they 
are spatial variations, these oscillations have a scale 
length of about 10 km. The spectral index n = -1.66 for the 
second spectrum in Figure 4 is nearly identical to one of 
the spectra in Figure 2. This near-equality suggests that 
similar transfer processes in spectral energy may be involved. 
(We note, however, that a specific spectral index does not 
uniquely define an instability process.) The third spectrum 
in Figure 4 corresponds to the. region of large peak-to-peak 
variations. The spectrum clearly shows the predominance of 
structure in the range 0.1-1.0 Hz (40-4 km). The spectral 
index is near -2, possibly due to the presence of steep 
edges in the irregularities. 

CONCLUDING COMMENTS 


The aurora-related F-region density irregularities 
discussed in this report illustrate the variety of plasma 
kinetic processes that are involved in high latitude ionospheric 
effects. Our study of the density irregularities is in a 
preliminary stage, and we expect that new and significant 
results will be forthcoming as the study advances. It is 
already obvious from an initial survey of available satellite 
data that aurora-related effects are detected outside the 
average auroral boundaries, i.e., both at higher and lower 
latitudes. Since the high latitude region is directly 
associated with magnetosphere-ionosphere coupling processes, 
we anticipate that the analysis of the ionospheric density 
irregularities in regions such as the polar cusp will address 
problems of global importance. 
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Ionosphere Equatorial irregularities 

"In Situ” Satellite measurements 


The S3-4 “in situ” satellite observations near the nighttime equator (2230 LT) have 
shown irregularity structures with depletions in plasma density up to 3 orders of magni¬ 
tude and ranges in horizontal extent from less than 1 km to tens of kms. The holes have 
sharper gradients on the eastern boundary than across the western counterpart, and the 
power spectral indices in the irregularities lie between 1.9-2.5 in the intermediate 

(Continues) 



















wavelength domain. The east-west assymmetry of the depletions is also shown in the 
irregularity intensity and spectral strengths. The variations in density are more intense 
on the western wall of the holes suggesting a scintillation enhancement on the western 
boundary. In addition, the results of power spectral analyses support the role of the 
Rayleigh-Taylor instability in the generation of intermediate scale size irregularities 
during the occurrence of equatorial spread-F. 
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TI1K STP/S3-4 SATELLITE EXPERIMENT: 
EQUATORIAL F-RKCION IRREGULARITIES 


INTRODUCTION 

Equatorial F-region irregularities have drawn much interest 
in recent years because of their undesirable effects on trans- 
ionospheric communications and their cause-effect relationship 
with fundamental plasma instability processes. In efforts to 
understand the causative mechanism(s), considerable advances have 
been made in areas of detailed ground-based radar observations 
and "in situ" measurements. The radar observations (e.g., 

Woodman and La Hoz, 1976; Tsunoda et al. 1979) have found that 
meter-size irregularities primarily populate the bottomside F- 
region in the early evening, at later times tend to rise up and 
break away from their lower altitude source regime and develop 
structures extending up to the 700-1000 km region. These structures 
have come to be called "plumes". 

On the macroscale (100's of meters to 10's of km) "in situ" 
nighttime equatorial measurements have revealed large biteouts in 
plasma density ranging up to three orders of magnitude (Hanson 
and Sanatani, 1973) and considered just as characteristic of 
spread-F as the much less intense (meter size) irregularities 
observed by radar. Later works on biteouts (Brinton et al. 1975; 
McClure et al. 1977) showed that ion composition inside and 
outside the holes can be vastly different. The molecular ions 
can be more abundant inside the holes than outside the holes, and 
the holes can vary from a few km to tens of km in the horizontal 
extent. 

In examining potential relationships between radar plumes 
and ionospheric depletions, Szuszczewicz (1978) suggested that 
equatorial holes and spread-F were the same phenomena with small 
scale irregularities imbedded within the large scale depletions. 

He argued that a chemical volume of ion density on the bottomside 
(containing the signature of bottomside species) could move 
upward through a stationary neutral atmosphere and appear as 
biteouts at higher altitudes with much smaller structures (down 
to the meter range) populating the density gradients which 
bounded the macroscale depletion. This model was in concert with 
the numerical results of Scannapieco and Ossakow (1976) and the 
drift measurements of McClure et al. (1977). 

To further study and definitively unfold the detailed 
relationships between large scale depletions, meter size 
irregularities and chemical transport processes, a coordinated 
investigation was conducted which involved simultaneous obser¬ 
vations by radar and "in situ" rocket-borne diagnostics 
(Szuszczewicz et al., 1980). The combined observations have 
shown that 
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(a) During conditions of well-developed equatorial spread-F 
the most intense "in situ" irregularities occurred on the bottom- 
side F-layer gradient. 

(b) Within a large scale topside F-layer depletion 
radar backscatter and "in situ" irregularity strengths 
maximized near the depletion's upper wall. 

(c) Ion composition within a topside depletion provided 
signatures of its bottomside source domain and estimates of 
average maximum vertical drift velocity. For long-lived 
depletions, it was found that molecular-ion signatures (NO 
and of) can be los£ while bottomside levels of N can be 
maintained when [0 3 £ N e >> [NO ] + [0 2 ]; and finally, 

(d) Large scale fluctuation^ of 0 + accompanied by a 
near-constant level of NO and of on the bottomside F-layer 
gradient suggests that neutral atmospheric turbulence was 
not a major source for bottomside ionospheric plasma irregu¬ 
larities and the associated triggering of equatorial spread- 
F. 


To complement the vertical profile information provided 
by the rocket observations, we present and analyze a sample 
of "in situ" measurements conducted on the STP/S3-4 satellite 
carrying a pulsed plasma probe experiment. The probe experi¬ 
ment employed a self consistent test for measurement integrity, 
while determining electron density and temperature as well 
as density and mean ion mass fluctuations at 5-20 meter 
resolution. The S3-4 experiment has been discussed by 
Szuszczewicz et al. (1981). In this paper, we discuss some 
of the recent results with emphasis on general horizontal 
morphology, relationships to basic instability processes 
and associated scintillation effects. 

EXPERIMENTAL RESULTS 
Plasma Depletions 

The data reported here was made available by the NRL-747 
paired-pulsed-plasma-probes (P4) experiment (Szuszczewicz et 
al. 1981) on STP satellite S3-4 in a sun-synchronous orbit 
at lower F-region altitudes. The satellite crossed the 
nighttime equator at about 2230 LT when the the frequency of 
occurrence of spread-F was high. One probe was biased to 
respond to variations in plasma electron saturation currents 
(I « N ) and the other probe responded to the ion saturation 








currents (I. = N /yfvL). Subject to the selection of a 
number of comman§able 1 modes of operation, either probe could 
be repetitively pulsed from its fixed-bias level using a 
special electronic procedure (Holmes and Szuszczewicz, 1975) 
to generate conventional Langmuir characteristics for full 
determinations of electron density N , temperature T and 
plasma potential V . The different modes of operation and 
experimental details are discussed by Szuszczewicz et al. 

(1981) . 

Figures 1 and 2 present samples of nighttime equatorial 
irregularity structures as measured by the currents collected 
by the electron (I ) and ion (I.) probes. The Figures show 
that the holes extlnd from a few kms to tens of kms, with 
depletion levels ranging from a factor of 3 (hole A) to a 
factor of 500 (hole F). A cursory analysis of Figures 1 and 
2 reveals three particularly interesting observations: 

(a) The density gradients on the opposing sides of each 
hole are different. Noting that the inclination of the 
satellite orbit is 96.4° (traveling in the east-to-west 
direction with time increasing left-to-right) it can be seen 
that density gradients in holes C-E (Fiq. 1) and F (Fiq. 2) 
are sharp on the eastern boundary and soft on each western 
counterpart. In the case of hole A, B and H (Fig. 1), the 
density gradients are not quite as sharp in comparison with 
holes C-F but the average density gradients are still softer 
on the western boundary. For the holes L (Fig. 1) and N' 
and L' (Fig. 2), the difference in density gradients on the 
opposing sides does not constitute any specific behavior 
while in the hole C' (Fig. 2) the density gradient is sharp 

on the western boundary. These observations can be summarized 
by noting that the density gradients are different across 
the boundaries of each depletion, with a preference for the 
density gradients to be sharper on the eastern side. 

(b) Another interesting feature of the holes is their 
similarities in structural morphology. Holes A, B and H are 
similar in their horizontal profiles. A similarity feature 
can also be identified in the depletion L (Figure 1) and L' 
(Figure 2). Both of these depletions have a center point (M 
and M' respectively) around which the structures look similar. 
At the central point (M and M'), the plasma density approaches 
the background (undisturbed) level suggesting that one large 
hole is breaking up or alternatively two smaller ones are 
merging. 

(c) The top panel of Figures 1 and 2 displays relative 
density fluctuations, <$I/I determined by variations about 


3 






PROBE CURRENT (AMP) 












Figure 2 

Horizontal profile of ionospheric F-region plasma 
density as indicated by electron (I ) and ion (I.) 
saturation currents measured on rev e #2470. 6I/I 1 

is the relative irregularity intensity as calculated 
by variations about linear detrends over sliding 
2.1 km intervals. 









linear detrends over sliding 2.1 km intervals throughout the 
depletions. The fluctuations (6I/I) easily identifies the 
degree of disturbance. The percent variation in the holes 
A, B and H (Figure 1) is less than 20% in most parts, while 
in other holes of Figures 1 and 2 the variation can be 
more than 50%. In addition, the variations can be different 
across the opposing boundaries of the holes. We discuss 
this more fully below. 


East-West Asymmetry 

Further illustrations of F-region irregularities appear 
in data revs 2122 and 2123 shown in Figure 3. Rev 2123 
shows four depletions (numbered 1 through 4) with the depletion 
level extending to more than two orders of magnitude. The 
rev 2122 shows three depletions (numbered 5 through 7) with 
depletions in density up to one order of magnitude. 

The discussion is facilitated by identifying certain 
features in rev 2123. First, there are clearly defined 
regions of undisturbed background ionosphere, marked alphabetically 
A through D; the smoothness of the relative density and the 
corresponding 0% fluctuations attest to their undisturbed 
nature. Focussing on depletions 3 and 4, we see that the 
irregularity intensities are 2 to 3 times larger on the 
western boundary than on eastern counterparts. This same 
relationship is true in depletions 1 and 2, but only after a 
qualification that suggests that 1 and 2 are halves of a 
larger depletion bounded by A and B. This is supported in 
part by the non-existence of a quiescient ionosphere between 
the two. When viewed from this perspective the western 
boundary is approximately twice as intense in irregularity 
intensity as the eastern boundary. In rev 2122, the depletions 
6 and 7 show the irregular intensity as more intense on the 
western boundary by a factor of 3-4. (On the basis of 
similar considerations applied to Figures 1 and 2, we note 
that the irregularity intensity is greater on the western 
boundaries of depletions D, E and H (Figure 1) and F (Figure 
2), while the eastern boundary in holes A-C (Figure 1) is 
more intense. Based on irregularity intensity (6I/I) 
observations, we conclude that there is a high probability 
of occurrence of more intense fluctuations across the western 
boundary of bottomside F-region depletions. 

The asymmetry in irregularity strengths and relationships 
to plasma instability mechanisms can be explored further 
through power spectral density analyses. We present in 
Figures 4 and 5 just such results for each of the boundaries 
in Figure 3 (IE and IW refer respectively to the eastern and 
western boundary of depletion number 1). Though the experiment 
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provides the density fluctuation power spectra with a maximum 
Nyquist frequency of 400 Hz (19 m resolution along the orbit 
and 3-5 m resolution perpendicular to magnetic field) the 
data in this figure has been decimated by a factor of 3 in 
order to include greater spatial domains across the depletion 
walls. Power spectral analyses are presented across the 
boundaries of each of the seven depletions with spectral 
indices (n, in the equation P = P f n ) ranging from 1.9 co 
2.5. More important however, is °the ratio of spectral strengths 
(entered as Pi hP q ( west)/P (east) in Figures 4 and 5) found to 
be 1.4 to 11.3 times more intense on the western boundaries. 

These intensity ratios extend down to a 15 meter wave¬ 
length perpendicular to the geomagnetic field. The spectral 
indices support the work of Keskinen et al. (1980) which 
predicts the same approximate range of values for horizontal 
irregularity structures perpendicular to B. The east- 
west asymmetry in the depletions is apparent from the 
linear detrend of data and the ratio of spectral strengths. 

We now show that the asymmetry bears on scintillation 
observations. To do this we note that scintillations depend 
upon AN(AN«AI) rather than AN/N (as calculated by 

(AI) and the power spectral analysis P Q « -(yi) 2 ). For the 
purpose of scintillation-effect calculations, we have determined 
|Ai| r -m.s. in fche wavelength domain 80 m to 8 km by assuming 

that P r ’ m ' s ‘ i n this domain is equal to (|Al| r * m " s ’/X) 2 . 

Defining P 2 as 

p _ (AI) r m s (in the wavelength range 80m-8km on the west wall) 

(AI) r m>s (in the wavelength range 80m-8km on the east wall) 

we find the values of P 2 ranging from 1.1 to 7.7. These 
results suggest enhanced scintillation effects on the western 
wall of the depletions, in agreement with the observations 
of Livingston et al. (1980). Furthermore, the radar measure¬ 
ments (Tsunoda, 1979) also show that the bottomside back- 
scatter strength is often asymmetric in the east-west plane with 
stronger backscatter from the western wall of a plume. The 
combined observations support a model of E-W asymmetry which 
allows for a neutral-wind driven instability growth rate 
enhancement on the western side of a rising bottomside F- 
region depletion. The maximum growth rate occurs on the 
western wall of the rising depletion, where the electron 
density gradient is most closely aligned with the plasma 
drift velocity vector in the reference frame of the neutral 
wind (Tsunoda, 1979; Zalesak et al. 1980). 





Figure U 


Power spectral analyses across the boundaries of the 
four depletions (1 through 4) of rev #2123 in Figure 3. 
The number in the upper right portion of each panel 
is the spectral index (n in the equation 

P = P Q k n ). The value P 1 in each right panel is the 
ratio of spectral strengths (P q ) of western to eastern 
boundary P = P (West)/P (East). The value P_ is the 
ratio of r.m.s. value of Spectral strength in the 
wavelength domain (80 m-8 Km) of western to eastern 

boundary (P 2 = p r , m .s. (West)/P r.m.s. (East)} ’ 
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Power spectral analyses across the boundaries of the 
four depletions (5 through 7) of rev #2122 in Figure 3 
The number in the upper right portion of each panel 
is the spectral index (n in the equation 

P = P k n ). The value P. in each right panel is the 
ratio 0 of spectral strengths (P q ) of western to easteri 

boundary P^^ = P (West)/P (East) . The value P_ is the 
ratio of rim.s. value o? spectral strength in the 
wavelength domain (80 m-8 Km) of western to eastern 

boundary (P, = P (West)/P (East)). 




COMMENTS AND CONCLUSIONS 


The results of the high resolution Si-4 satellite 
experiment for equatorial F-region irregularities show large 
scale plasma depletions (1 km to 10*s of km wide) with 
smaller scale irregularities superimposed {smallest detectable 
scale size ^ 10 meters perpendicular to B). The depletions 
(which may in fact represent the bottomside upwelling process 
that has come to be identified with the lower F-region 
manifestations of spread-F) show east-west asymmetry with 
the irregularity intensity and spectral strengths generally 
more intense on the western boundary. Associated calculations 
over density fluctuations in the range 80 meters to 8 km 
suggest that scintillation effects would be similarly more 
intense on the western boundary. We find these observations 
consistent with radar (Tsunoda, 1979) and scintillation 
measurements (Livingston et al. 1980) as well as the recent 
computational work of Zalesak et al. (1980). Furthermore, 
we find the horizontal power spectral indices to lie between 
1.9-2.5, supporting the role of Rayliegh-Taylor instability 
(Keskinen et al. 1980) and the recent theoretical and experi¬ 
mental comparison in the vertical plane (Keskinen et al. 

1981) . 
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Ionospheric irregularities 
Power spectrum 
Satellite measurements 

Plasma diagnostics | 

The pulsed plasma probe technique has been expanded to include simultaneous determinations of 
absolute electron density, density fluctuations, electron temperature, and mean ion mass with resolution 
limited only by probe geometry, sheath size, and telemetry. The technique has been designed to test for 
coupling of electron density variations and ion composition irregularities in multi-component plasmas 
by the comparison of electron density fluctuation power spectra Pjq(k) and a newly developed 
diagnostic parameter, the mean ion mass fluctuation spectra -* P M (k). In addition, the experiment 
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extends satellite-borne irregularity spectral analyses down to the 5-20 meter range while attempting to 
synthesize F-region plasma instability processes on the basis of N e , T e , AN e , P N and P M . Initial 
results demonstrate the expanded diagnostic capability for high spatial resolution measurements of 
mean ion mass and provide experimental evidence for the role of ion composition in multi-stepped 
plasma instability processes. Specific results include a spectral index X n in P N - A n f* x n of 1.6-2.9 
over the wavelength range from 1 km to 6 meters under conditions identified with an unstable equatorial 
nighttime ionosphere. Simultaneous measurements of SMj/Mj •* (Pj^ - A m f x m) and 5N e /N e -* (P N 
- Anf*^) have shown a general behavior tending to lower power (A m < A n ) and softer spectra (X^ 

< X,,) in ion mass fluctuations when compared with fluctuations in total plasma density. Limited 
analyses of the two power spectral elements raises hopes for the differentiation between plasma 
mechanisms that can lead to similar indices in p N- 
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THE S3-4 IONOSPHERIC IRREGULARITIES SATELLITE EXPERIMENT: 
PROBE DETECTION OF MULTI-ION COMPONENT PLASMAS 
AND ASSOCIATED EFFECTS ON INSTABILITY PROCESSES 


In recent years interest 


ionospheric research ha 


turned from the definition of the zero-order ionosphere to 
an understanding of the more dynamic processes involving 
plasma instability mechanisms and their cause-effect relation¬ 
ships. While much progress has been made, more detailed 
experimental and theoretical work is needed to unfold the 






























































(Szuszczewicz et al., 1980; Narcisi and Szuszczewlcz, 1981), 
and cascade the irregularity distributions from 10's of 
kilometers to fractions of a meter (e.g., Keskinen et al.,1981 


Kelly et al., 1981). 

Outside the latitudinal domain of equatorial spread-F, 
the regions of primary interest in the S3-4 study of F- 
region irregularities involve the main trough, auroral oval, 
and domains encompassing the ring current, polar wind and 
the cusp. The mid-latitude and dayside equatorial F-regions 
are generally very regular in structure and consequently of 
less interest. 

Within the framework of the relatively simple picture 
in Figure 1, the improved plasma diagnostic capabilities on 
the S3-4 satellite and associated analysis efforts will 
attempt to assemble a relatively comprehensive catalogue of 
polar and equatorial irregularities along with an improved 
understanding of their cause-effect relationships. Indeed, 
some initial efforts have already been published (Singh et 
al., 1981; Rodriguez et. al., 1981; and Szuszczewicz et al., 
1981). In the subsequent sections, the experimental technique 
will be described, associated data presented, and evidence 
will be provided which suggests that ion composition plays 
an important role in multi-stepped plasma instability 


processes active in the F-region ionosphere. 
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Fig. 1 — Schematic representation of F-region plasma density irregularities and associated geophysical domains 










THE EXPERIMENTAL TECHNIQUE 


The NRL-747 experiment employs a pair of pulsed plasma 
probes (P^ is the designated acronym), each of which is 
capable of simultaneous measurement of electron density, 
temperature and density fluctuation power spectra regardless 
of the state of turbulence or the degree of irregularity in 
the ionospheric plasma medium. Together, the pair of probes 
also provide mean ion mass fluctuation measurements to a maximum 
Nyquist frequency of 200 Hz. The instrument is a Langmuir- 
type probe using a special electronic procedure for generating 
the current-voltage characteristic (Holmes and Szuszczewicz, 
1975; Szuszczewicz and Holmes, 1975). The pulse modulated 
approach reduces the distortion of the current-voltage 
characteristic that can result from surface contamination 
and allows millisecond tracking of density fluctuations that 
might occur during the time required to generate the current 
voltage characteristic. (The temporal resolution is limited 
only by telemetry.) 

Subject to the selection of 1 of 8 commandable modes of 
operation, each of the probes had applied to it some variation 
of the voltage function illustrated in Figure 2. The pulse 
modulated waveform, following the sawtooth envelope, provides 
the fundamental data set for a "conventional" Langmuir 
current-voltage characteristic and associated determinations 
of N g and T g (Chen, 1965) at a nominal 3 Hz rate. During 
the interpulse period, a fixed-voltage is applied, with 
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population, T 


is the associated temperature of an assumed 


e (i ) 

Maxwellian distribution, e is the fundamental electron 

charge, k is Boltzmann’s constant, w is the satellite velocity 
e (i) 

and <f p is the baseline voltage Vg applied to the E (I) 

A ( 4 \ 

probe and referenced to the plasma potential ($^ ■ Vg- V B ). 

The square of the ratio I £ (sat)/I i (sat) can be written 
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For laboratory and rocket-borne experiments equation (4a) 
would apply, whereas in the S3-4 satellite investigation 
equation (4c) applies. Equation (4c) is plotted in Figure 
3 for two sets of bias potentials, (|4> a |, l4*p|) “ (2V,1V) 

and (IV,2V). The results in Figure 3 show that with limitations 
placed on the ratio ^p/^p* variations in (l^/I^) 2 can be 
taken as a direct measure of ion mass variations. 












DATA ANALYSIS AND RESULTS 


Initial reduction routines for bulk processing and 
plotting of P^/S3-4 data begin with an orbit-by-orbit plot 
of relative electron density as measured by changes in ion- 
and electron-saturation currents. A representative sample 
of this data collected on orbit 2177 is shown in Figure 4 
where the abscissa coordinates are universal time (UT), 
altitude (ALT in km), latitude (LAT), longitude (LONG), 
magnetic latitude (MLAT), and L-shell value (L). The probes' 
magnetic aspect angle is also plotted in the figure. 

The left hand edge in Figure 4 corresponds to the 
satellite's ascending node (south-to-north) in the midnight 
hemisphere near the south magnetic pole (MLAT ■ -90°). 

With increasing time (UT) the satellite passes through the 
nighttime equator, the main trough, over the northern auroral 
oval and into the dayside ionosphere where vehicle solar cell 
voltage biases the entire vehicle such that both probes draw 
approximately equal ion saturation currents. 
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It is appropriate to note that the experimental determination 

of mean-ion-mass fluctuations SM^(-*-P^), through variations 
o 

in fl e (sat)/I i (sat)] , assumes the relative constancy of all 

potentials. (This includes the spacecraft potential as well 
as the potentials which each probe presents to the plasma.) 

The spacecraft potential can vary as a direct result of 
changes in local plasma density since the floating potential 
of a body is dependent upon the ratio of its radius to the 
Debye length. For large space vehicles however, floating 
potential variations caused by even substantial plasma 
density variations should be small (Szuszczewicz, 1972). 

Another possible source of potential variations involve 
charging of contamination layers on the vehicle and/or on 
the probes (Szuszczewicz and Holmes, 1975). From the S3-4 
data examined to date, variations in (I e /I^)^ associated 
with charging on contamination layers, appear to be a slowly 
varying function of time with no attendant effects on P^. 
Therefore, we have concluded that the spectral dependence of 
P M is indeed representative of. variations in mean-lon-mass 
48 i- 

To experimentally demonstrate P^(k) and the associated 
relationship with P N (k) consider the high-resolution measure¬ 
ments (rev #2123) of the relative electron density across 
the nighttime equator (Figure 5). The peak electron density 
is approximately 5 (10^) cm ^ at I g - 3 (10 ®) amp. The 
large scale depletions are seen to extend to 2 orders of 
magnitude with widths ranging from 50 to 170 km over a 600 km 
orbital segment. 
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TIME SECONDS 


UT 10:50:22 
ALT 241.1 km 
LAT-18.71 °S 
LONG 175.9°E 


10:58:42 
219 km 
11.9°N 
170.3°E 


Fig. 5 — An expanded view of relative electron density encountered during the nighttime equatorial 
crossing on 83-4 rev 2123. The relative electron density is presented by baseline electron saturation 
current. 
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6Mi/Mi is a fairly complicated function of M^/Mg, 

N°/N° N 1 /No and 6N /N itself (See Eq. 7b); and at this 
a p ot p e e 

point we can only speculate on the many manifestations that 
and P M might take for the varied ionospheric conditions 
encountered in the S3-4 mission. For example, it has been 
suggested (J.R. Goldman, private communication) that differences 
in gradient scale lengths for and g's would result in 
a more rapid fall off with increasing k for the quantity 
with an initially larger gradient scale length. This difference 
should be a direct observable through the P N and P M determination. 
Furthermore, there is the possibility that the simultaneous 

_2 

measurement of P M (k) could help differentiate between a k 

_2 

spectrum due to sharp edges and a k spectrum due to 
gradient-drift or say drift-dissipative waves. 

Currently, we have not analyzed sufficient data to 
exact any of the above analyses. We are at the beginning of 
a synoptic study designed to systematically unfold the 
interdependence of P N and P^, and the contributing roles of 
ion mass distributions through M^/Mg and N°/Ng. However, 
the results in Figure 6 as well as in Figure 7 can be used 
to reflect some of our observations: 

(i) The absolute power in P N is generally much greater 
than that in P^. In Figure 6, for example, power law fits 
yield A Q / A m ^ 86, reflecting a situation in which the 
absolute density fluctuations 6N g /N e are two orders of 
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magnitude more Intense than counterparts in 

(ii) More often than not the k-dependence of P N is 

maintained to shorter wavelength (higher frequencies) than P M (k). 

(ill) The relationship between.^the spectral Indices of 

P N and P M is highly variable. Figure 7 is an illustration 

of a random selection of 1.2-second interval spectra near 

the domain sampled in Figure 6. Panel A represents the power 

spectrum of the quiescent region before the depletion, while 

panels B through D represent contiguous power spectra within 

the depletion. In the quiescent region (Panel A), both 

density and mean ion mass fluctuations show noise spectra, 

while inside the depletion well-defined k-dependencies in 

P N and (less frequent) k-dependent spectra in P M are observed. 

While small in number, this sample of data is rather representative 

of results accumulated to date. More frequently than not P M 

is equal to or softer than P N (e.g., panels B,C and D). In 

addition, there are observations of well-defined k-spectra 

in P N with none in P M (i.e., only noise, as in panel D). 

Panel D Is interpreted as an unstable plasma environment 

with a single component ion population. This is In agreement 

with the predictions of Eq. (7b) with M - M„. 

a p 

Results like those outlined above are continuing with a 
view to a synoptic definition of the PP N interdependence 
and the associated roles of M^, gradient scale lengths and 
geoplasma conditions (N g ,T e and MLAT). 


22 






f th 














experiment goes beyond the 70 meter limit of earlier satellite 
measurements and lends itself to additional arguments which 
support the association of the Rayleigh-Taylor instability 
occurrence of equatorial spread-F. 

(2) Measurements of 6M^ with resolution that allows 
detection of macroscale features (e.g., transitions from 0 + 
dominant to N0 + dominant environment) and microscale spectral 
analysis down, to dimensions comparable to the 6N g capability. 
As expected, power law fits (P (k) - Af X ) to P N and P^ have 
revealed considerable variations with X M generally (but not 
always) smaller than X N . Similarly, the absolute power A M 
in ion mass fluctuations have been found to be considerably 
less (@ 2-100) than the power in P N . A systematic study of 
the detailed relationships and dependence on zero-order 
plasma conditions is currently underway. 
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